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ARE GIANT PLANETS FORMING AROUND HR 4796A? 
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ABSTRACT 

We have obtained FUSE and HST STIS spectra of HR 4796A, a nearby 8 Myr old main sequence star 
that possesses a dusty circumstellar disk whose inclination has been constrained from high resolution 
near-infrared observations to be — 17° from edge-on. We searched for circumstellar absorption in the 
ground states of C II A1036.3, O I A1039.2, Zn II A2026.1, Lyman series H 2 , and CO (A-X) and failed to 
detect any of these species. We place upper limits on the column densities and infer upper limits on the 
gas masses assuming that the gas is in hydrostatic equilibrium, is well-mixed, and has a temperature, 
Tgas ~ 65 K. Our measurements suggest that this system possesses very little molecular gas. Therefore, 
we infer an upper limit for the gas:dust ratio (<4.0) assuming that the gas is atomic. We measure less 
gas in this system than is required to form the envelope of Jupiter. 

Subject headings: stars: individual (HR 4796A) — circumstellar matter — planetary systems: formation 



1. INTRODUCTION 

Planets, asteroids, and comets are believed to form in 
circumstellar disks around stars with ages <few 100 Myr. 
Current (uncertain) models suggest that Jovian planets 
form either via rapid gravitational collapse through disk 
instability within a few hundred years (Boss 2000) or via 
coagulation of dust into solid cores within the first — 1 
Myr and accretion of gas into thick hydrogen atmospheres 
within the first -10 Myr (Ruden 1999). At present, the 
timcscales on which giant planets form and accrete their 
atmospheres have not been well constrained observation- 
ally. Studies of objects with ages —10 Myr afford the op- 
portunity to constrain the timescale for gas dissipation 
around young stars and therefore the timescale for giant 
planet formation. 

Although the circumstellar disks around pre-main se- 
quence stars possess both gas and dust, the circumstellar 
disks around main sequence stars may be gas poor. Mil- 
limeter searches for CO around stars with ages 1-10 Myr 
suggest that within a few Myr, less than a Jupiter mass 
of molecular gas remains in the circumstellar environment 
(Zuckerman, Forveille, & Kastner 1995). However, CO 
may be depleted in the outer regions of disks because it 
freezes out onto dust grains at temperatures, T gas ~ 20 
K. Since gas with interstellar or solar abundance is mainly 
composed of hydrogen and H 2 does not freeze onto grains, 
recent searches for circumstellar gas have focused on H 2 . 
ISO searches for emission from the S(0) and S(l) tran- 
sitions of H 2 at 28 (j,m and 17 /jm respectively, suggest 
that 50 - 2000 M ffi H 2 exists around j3 Pic, 49 Cet, and 
HD 135344 (Thi et al. 2001). However, FUSE searches 
for H 2 absorption in the Lyman band toward (3 Pic place 
upper limits on the column density of H 2 , N(H 2 ) <10 18 
cm -2 , which is 3 orders or magnitude lower than that in- 
ferred from ISO (Lecavelier des Etangs et al. 2001). The 
Thi et al. (2001) results not only apparently conflict with 
the ultraviolet H 2 measurements but also with millimeter 
emission line studies which suggest that j3 Pic possesses 
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<10~ 5 M e CO (Liseau & Artymowicz 1998). Ground- 
based follow-up of ISO H 2 detections toward Herbig Ae 
and T-Tauri stars yielded upper limits which were smaller 
than ISO detections (Richter et al. 2002), suggesting that 
the H 2 must lie in an extended cloud. However, an H 2 
cloud around (3 Pic should have produced a detectable ul- 
traviolet absorption line (Jura 2003). Recently, fluorescent 
H 2 Lyman-band emission has been detected from the sur- 
face layers of the circumstellar disk around the —10 Myr 
old star TW Hydrae (Herczeg et al. 2002) using HST 
STIS. We have carried out a FUSE and STIS search for 
circumstellar molecular and atomic gases around the —8 
Myr (Stauffer et al. 1995) debris disk object HR 4796A, 
another member of the TW Hydrae Association, to help 
constrain the gas dissipation timescale. 

HR 4796A is a main sequence A0V star with a Hip- 
parcos distance of 67.1 pc (see Table 1 for a summary of 
stellar properties) and a fractional dust luminosity L/^/L* 
= 5xl0~ 3 (Jura et al. 1998). Submillimeter observations 
suggest that it possesses >0.25 M e dust (Greaves et al. 
2000), more circumstellar dust than any other main se- 
quence A-type star in the Yale Bright Star Catalog (Jura 
et al. 1991). High resolution scattered light imaging has 
resolved a narrow, dusty ring at 70 AU from HR 4796A 
with an inclination of —17° from edge-on (Schneider et 
al. 1999). High resolution mid- infrared imaging has re- 
solved the disk in thermal emission (Koerner et al. 1998; 
Jayawardhana et al. 1998). The presence of a central 
clearing (Jura et al. 1998) and of a brightness asymmetry 
from the two lobes of the disk (Telesco et al. 2000; Wein- 
berger, Becklin, & Schneider 2000) have led to speculation 
that the circumstellar disk possesses a planetary compan- 
ion at <70 AU (Wyatt et al. 1999). Blackbody fitting 
to the IRAS fluxes suggests that the grain temperature, 
T d ust = 110 K (Jura et al. 1995). 

Previous studies of circumstellar gas around HR 4796A 
have not been very successful. Searches for millimeter 
and submillimeter emission from CO (J=3-^2, J=2— 
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J=l — >0) , used as a tracer for H 2 , have failed (Zuckerman 
et al. 1995, Liseau 1999, & Greaves et al. 2000), per- 
haps because the angular size of the source as measured 
in the infrared is significantly smaller than the beam of 
available telescopes. Another possibility is that circum- 
stellar molecular gas is not present because it is rapidly 
photodissociated under the intense stellar ultraviolet field 
of the central AOV star. Holweger, Hcmpel, & Kamp 
(1999) have detected two narrow, circumstellar Ca II K 
absorption features, one component close to the stellar ve- 
locity and another weaker one blueshifted by ^10 km/sec 
which they attribute to possible circumstellar gas. 

2. OBSERVATIONS 

HR 4796A was observed with the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) for 12.4 kscc on 2002 April 25 
in histogram mode, using the low resolution 30" x 30" 
aperture which covers the wavelength range between 905 
and 1187 A. This wavelength coverage allows us to search 
for circumstellar absorption in the ground states of C II 
A1036.3 and O I A1039.2, and in Lyman scries H 2 . The 
FUSE satellite has four channels (LIF 1, SIC 1, LIF 2, 
SIC 2) which form two nearly identical "sides" (labeled 1 
and 2) that each consist of a LIF grating, a SIC grating, 
and a detector. Each detector is divided into two indepen- 
dent segments (A and B), seperated by a small gap. The 
eight partially overlapping spectra that fall on different 
portions of the two detectors (LIF 1A, LIF IB, etc.) cover 
the entire wavelength range. A description of the on orbit 
performance of FUSE is described in Sahnow et al. (2000) 
with a conservative estimate for the spectral resolution (R 
= A/AA = 15000) of the spectrograph. The data were 
calibrated at Johns Hopkins using the CALFUSE 1.8.7 
pipeline. 

HR 4796A was observed with Space Telescope Imaging 
Spectrograph (STIS) in echelle mode on the Hubble Space 
Telescope (HST) for 930 sec and 550 sec on 2002 August 
22, with the E140H and E230H gratings respectively using 
the 0.2"x0.09"slit. The spectra covered the wavelength 
regions 1425 A - 1620 A and 1880 A - 2155 A with a spec- 
tral resolution (R = A/AA = 220,000). This wavelength 
coverage allows us to search for circumstellar absorption 
in the ground state of Zn II A2026.1 and in CO (A-X). 
The data were calibrated at the Space Telescope Science 
Institute using the standard STIS software package calstis, 
executed under the OPUS pipeline. 

We observe absorption in the ground states of Si II at 
1526.14 A, Fe II at 1608.45 A, and Zn II at 2026.14 A along 
the line of sight to HR 4796 A (see Figure 1). The Si II 
and Fe II lines possess two absorption components, one at 
about -5 km/sec and another at about -15 km/sec, con- 
sistent with the Ca II results from Holweger et al. (1999). 
We believe that both gas components are interstellar. Cir- 
cumstellar gas in Keplerian orbit at 70 AU from an AOV 
star would have a velocity of ~5.5 km/sec. Therefore, any 
circumstellar gas features will be indistinguishable from 
interstellar lines on the basis of line width. However, cir- 
cumstellar gas at 70 AU from an AOV should be radiatively 
pumped, populating the excited fine structure lines. We 
observe no absorption in the Si II* line at 1533.43 A, sug- 
gesting that the Si II gas is interstellar. Since the Fe II line 
has the same velocity structure as the Si II line, this ma- 



terial is probably also interstellar. Even though the Zn II 
line only possesses one velocity component, its velocity is 
similar to that of the Si II and Fe II components near the 
stellar velocity. 

We observe no absorption in the H 2 Lyman (5,0) band 
or the CO A-X (1,0) band toward this star. The 3cr upper 
limits on the H 2 line-of-sight column densities in the J = 
and J = 1 rotational levels are measured, from the P(l) 
and R(0) transitions at 1036.54 A and 1038.16 A, to be 
N(J = 0) < 10 15 cm" 2 and N(J = 1) < 3.7 x 10 15 cur 2 
using wavelengths and oscillator strengths from Abgrall et 
al. (1993). The 3<r upper limits on the CO line-of-sight 
column densities in the J — 1 and 3 — 2 rotational levels arc 
measured, from the R(0) and R(l) transitions at 1544.45 A 
and 1544.39.73 A, to be N(J = 1) < 1.2 x 10 13 cm" 2 and 
N(J = 2) < 2.4 x 10 13 cm~ 2 using wavelengths and oscil- 
lator strengths from Morton & Noreau (1994). Similarly, 
we observe no absorption in C II A1036.3 and O I A1039.2, 
and no circumstellar absorption in Zn II A2026.1 a tracer 
for atomic hydrogen in the interstellar medium (York & 
Jura 1982). The 3<r upper limits on these column densi- 
ties are N(C II) < 1.4 x 10 14 cur 2 , N(0 I) < 2.8 x 10 15 
cm~ 2 , and A(Zn II) < 2.6 x 10 12 cm~ 2 and are estimated 
using oscillator strengths from Morton (1991). We calcu- 
late equivalent width upper limits assuming circumstellar 
line widths of 5.5 km/sec for the STIS data. We calculate 
equivalent width upper limits assuming circumstellar line 
widths of 2 resolution elements, after binning the data to 
a resolution of 0.05 A or ^29 km/sec for the FUSE data. 

3. DISK VERTICAL STRUCTURE 

Circumstellar gas in the low density environments 
around Vega-type stars is primarily heated via photoelec- 
trons ejected from dust grains by stellar ultraviolet pho- 
tons and primarily cooled by radiation from the [C II] fine 
structure line at 157.5 /im (Kamp & van Zadelhoff 2001). 
HR 4796A possesses an M-type companion, located 7.7" 
away from the star. The companion has a ROSAT 0.1 - 2.4 
keV luminosity, L x = 1.5 x 10 29 ergs/sec (Jura et al. 1998) 
which corresponds to an X-ray flux of F x < 2 x 10~ 4 ergs 
s _1 cm~ 2 at the position of HR 4796A. We estimate the 
photoelectric and X-ray heating rates for the gas around 
HR 4796A to determine whether X-ray heating plays an 
important role in this circumstellar environment. 

The photoelectric heating rate for micron-sized silicate 
grains, like the ones found around HR 4796 A (Augereau 
et al. 1999; Li & Lunine 2003), is 

T pe = 2.5 x 10- 4 ae X n H (1) 

(Kamp & van Zadelhoff 2001) where nn is the number 
density of hydrogen, a is the grain cross section per hy- 
drogen nucleus, e is the photoelectric efficiency, and \ is 
the 912 - 1110 A photon flux measured in units of the 
Habing field (F H = 1.2 x 10 7 cm~ 2 s _1 ). At a distance 
of 70 AU from HR 4796A, we estimate \ = 460 from our 
FUSE spectra, roughly consistent with a Kurucz stellar 
atmosphere with T e ff = 10,000 K, logg = 4.0, and solar 
metallicity (Fajardo-Acosta, Telesco, & Knacke 1998). For 
large silicate grains, the dust UV cross section a = 2.34 x 
10 -21 /<5 cm 2 H-atom^ 1 where 5 is the gas:dust ratio and 
e = 0.06 for low temperatures. If the gas:dust ratio is 4 (a 
= 5.85 x 10~ 22 cm 2 H-atouT 1 ), e = 0.06, and X = 460, 
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then the photoelectric heating rate is T p 



4.0 x 10 



-24 



tih ergs cm s . 

The X-ray heating rate for a cold neutral gas is 

= h>?xF x n H (2) 

where //,(~ 0.3) is the fraction of the absorbed X-ray en- 
ergy that heats the gas, a x is the X-ray cross section due 
to all atoms per hydrogen nucleus, and F x is the X-ray 
flux. If the bulk of the X-ray energy is carried by 0.5 keV 
photons, then the cross section a x = 1.65 x 10~ 21 cm 2 
(Maloney, Hollenbach, & Tielens 1996). If F x < 2 x 10~ 4 
ergs s _1 and a x = 1.65 x 10~ 21 cm 2 , then T x < 9.9 x 
10~ 26 nn ergs cm~ 3 s _1 , suggesting that X-ray heating is 
not significant in the optically thin regime. X-ray heating 
in more gaseous circumstellar disks is even less important 
because the disk becomes optically thick to the X-ray ra- 
diation. 

Kamp & van Zadclhoff have derived an analytic approx- 
imation for the gas temperature assuming that the gas is 
heated by photoelectrons ejected by dust grains and is 
cooled via [C II] line emission at 157.7 fim. In this case, 
the gas temperature 

T 9 -l^ (3) 



log[2( 



e c A 10 hv 10 
2.5xl0- 4 e<rx 
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where ec is the carbon abundance, A w and hvio are the 
Einstein A coefficient and energy carried by the photon 
in the C II transition. For the 157.7 /im [C II] line, Aio 
= 2.4 x 10~ 6 s- 1 and hv w = 1.27 x 10~ 14 ergs. If the 
circumstellar environment around HR 4796A possesses a 
solar carbon abundance, then ec = 4.6 x 10~ 4 and the 
gas temperature T ffas = 65 K. If the carbon abundance 
is interstellar ec = 1.4 x 10~ 4 , then the cooling via car- 
bon emission will be less efficient and the circumstellar gas 
will be hotter. In Figure 2, we plot the gas temperature 
as a function of the circumstellar gas mass (or gas:dust ra- 
tio) as inferred from equation (3). For HR 4796A, this 
approximation only holds for the outer portions of the 
disk. For example, O I is an important coolant at the 
inner edge of the disk midplane. In addition, this approx- 
imation breaks down for HR 4796A models with gas:dust 
ratio larger than a few. In higher gas:dust ratio models, 
the implied gas densities are large enough to allow CO to 
become self shielding, unlike the circumstellar disk models 
around Vega, because the HR 4796A circumstellar disk is 
more compact. In these cases, the cooling also becomes 
more complicated. 

We assume that the circumstellar gas is in hydrostatic 
equilibrium. In this model, the gas number density, n, 
can be written as a function of the height above the disk 
midplane, z, 

n(z) = n oe ^(t) 2 where| = (^|^)i (4) 

where H is the disk thickness parameter, R is the distance 
to the star, M* is the stellar mass, and fj, is the mean 
molecular weight of the gas. If the gas is primarily H 2 
with [He]/[H 2 ] = 0.2 so that fi = 3.9 x 10~ 24 g. If the 
gas is primarily H I with [He]/[H2] = 0.1, then fi = 2.2 
x 10~ 24 g. Since HR 4796A is not viewed edge-on, suc- 
cessful detection of circumstellar gas (or the establishment 
of meaningful gas mass upper limits) via absorption line 
studies requires that the disk thickness parameter, H ~ 



R sin i. Smaller thickness parameters lead to less puffed- up 
disks which provide fewer gas atoms and molecules along 
the line of sight. The inclination and the radius of the HR 
4796A disk has been measured to be ~17° from edge-on 
and 70 AU (Schneider et al. 1999) suggesting that a disk 
thickness parameter of ~10 AU is necessary. If the gas 
is molecular, M* = 2.4 M Q , and T gas — 65 K, then the 
thickness parameter H — 6.1 AU at R = 70 AU and our 
line of sight intercepts the gas at 3.4 H above the disk. If 
the gas is atomic, then the thickness parameter H — 8.1 
AU at R = 70 AU and our line of sight intercepts the gas 
at 2.5 H above the disk. 

4. GAS MASS ESTIMATES 

We estimate the gas mass upper limits from our line-of- 
sight column density upper limits for molecular and atomic 
gases assuming that the gas species are well-mixed and in 
hydrostatic equilibrium, (i.e. If the gas is molecular, we 
will assume a disk thickness parameter H = 6.1 AU for 
both CO and H 2 despite the fact that CO molecules are 
more massive than H2 molecules. If the gas is atomic, 
we will assume a disk thickness parameter H = 8.1 AU 
for C II, O I, and Zn II.) To estimate the circumstellar 
gas mass from the column density upper limits, we must 
assume a particular gas mass surface density, S. For sim- 
plicity, we assume that the gas and dust have the same 
radial distribution. If the dust and gas are not cospatial 
then photoelectric heating will not effectively warm the 
gas and the gas temperature will be significantly colder 
than 65 K. 
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(Klahr & Lin 2000) where R = 70 AU is the peak of the 
dust distribution, 5R = 10 AU, and d = 0.9 (which cor- 
responds to Ri = 68.1 AU). We can write an expression 
for the gas density from equations (4) and (5) 

E(r) 

n(r,z)= y > - e *W>> (6) 
V27rii(r) 

We estimate the column density by integrating the mass 
density in equation (6) along the line-of-sight. If the line 
of sight is inclined i° from the disk midplane, then the 
estimated column density, N, is 



N 
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g(g) 

m s cos i^/2n JR in H(R) 



e dR (7) 



where m s is the mass of the species of interest, i = 17°, 
T gas = 65 K, M* = 2.4 M Q , R in = 60 AU and R out = 80 
AU. 

If we set the estimated column density to be less than 
or equal to the observed column density upper limit, we 
can establish an upper limit for So. Once So is known, 
estimating the circumstellar gas mass upper limit only re- 
quires integrating the mass surface density. For molecular 
gas around HR 4796A, the gas mass in species s can be 
written as the following expression: 



M gas < 2.9 x 10" 
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where N is the measured column density. The J = and 
J = 1 H2 column density upper limits listed in Table 3 
suggest H 2 gas masses of <2.5xl0~ 4 M® and <9.1xlCT 4 
M®, respectively. This result suggests that the HR 4796 A 
circumstellar disk possesses very little molecular gas. 

The nondetection of H 2 around HR 4796A may not 
be surprising because HR 4796A, an AOV star, pro- 
duces a high ultraviolet photon flux which is capable of 
rapidly photodissociating molecular gas. Therefore, we 
also searched for absorption due to circumstellar O I, C II, 
and Zn II, a tracer for atomic hydrogen in the interstellar 
medium (York & Jura 1983). We can estimate the gas 
mass assuming that the bulk of the gas is predominately 
atomic. 

We list gas mass upper limits for each of the atomic species 
in Table 3 using equation (9). We also infer atomic hydro- 
gen disk mass upper limits assuming that the composition 
of the circumstellar gas is approximately solar. 



M gas < 3.3 x 10" 



M H < 1.4 x 10~ 8 - ( N _ 2 
e s V 10 i2 cm 1 



(10) 



where e s is the abundance of the species of interest. For 
e = 8.1 x 10~ 4 , e c = 4.6 x 10~ 4 , and e Zn = 3.8 x 10~ 8 
(Gray 1992), we find M H <1.0 M®. Since submillimeter 
continuum measurements estimate >0.25 M® dust around 
HR 4796A (Greaves et al. 2000), our hydrogen gas mass 
upper limit corresponds to a gas:dust mass ratio of <4.0. 

If circumstellar gas is in hydrostatic equilibrium, then 
any gas mass estimates are dependent on the inferred gas 
temperature. The density depends on the exponential of 
-(1/H) 2 ; therefore, small changes in the disk thickness pa- 
rameter can lead to dramatic changes in the inferred gas 
mass. We plot the inferred H I mass from our Zn II column 
density upper limit as a function of gas temperature in Fig- 
ure 3. For example, changing the gas temperature from 65 
K to 28 K corresponds to changing the disk thickness pa- 
rameter from 8.1 AU to 5.3 AU and the inferred H I mass 
from <1.0 M® to 40 M®. This corresponds to a gas:dust 
ratio <160, approximately the interstellar gas:dust ratio. 

5. THE GAS TEMPERATURE AND DISK CHEMISTRY 

Kamp & van Zadelhoff (2001) have produced detailed 
numerical models of circumstellar gas heating and cooling 
around main sequence A-type stars to calculate the gas 
temperature and chemistry. Given a prescribed gas and 
dust density distribution, these models solve the chem- 
istry and energy balance of the gas self-consistently on a 
two dimensional grid. The computation proceeds along a 
number of rays originating at the star's position under var- 
ious angles. Typically, 500 radial gridpoints and 30 angles 
covering about 3 scaleheights are considered. The heating 
in low density disks, like the one found around HR 4796A, 
is dominated by photoelectrons ejected from circumstellar 
dust grains by stellar ultraviolet photons. Kamp & van 
Zadelhoff (2001) use ATLAS9 models to estimate stellar 
photospheres and assume that circumstellar dust grains 
are large silicate grains with radius, a — 3 /im, and density, 
p = 3.0 g cm~ 3 . They assume a dust UV extinction cross 
section ajjv — 2.34 x 10~ 21 /<5 cm 2 (H-atom) -1 and a fixed 
fraction of vibrationally excited H 2 , fm = 1-0 xl0~ 5 . 



We have improved the Kamp & van Zadelhoff (2001) 
models by making the following changes: (1) including 
cosmic ray reactions, (2) lowering the temperature for CO 
freeze out to ^20 K, (3) including an escape probability 
formalism for line photons (Tielens & Hollcnbach 1985; 
TH85), and (4) calculating the statistical equilibrium for 
C II following the approach used for O I and CO by Kamp 
& van Zadelhoff (2001). Cosmic ray reactions are added to 
make the code more general. This change does not affect 
the result described here because the chemistry around 
main sequence A-type stars is driven by stellar photons. 
The escape probability formalism is necessary because the 
main cooling lines, such as [O I] and CO, become opti- 
cally thick for disks with gas:dust ratios of 10 and 100. 
The escape probability of a line photon, (3, is conserva- 
tively estimated from the line optical depth towards the 
star at each gridpoint. Since the optical depth perpendic- 
ular to the disk is much lower than the optical depth in the 
disk midplane, our estimates yield reasonable lower limits 
for the line cooling and therefore reasonable upper limits 
for the gas temperature. We also include a model which 
does not contain the one dimensional escape probability 
approximation to estimate the possible error arising from 
this approach. 

We have modeled the circumstellar gas around HR 
4796A using the model described above with a gas den- 
sity distribution, 



n{r) = Hi ( — 

. Hi 



-2.5 



(11) 



where Ri is the inner disk radius and rij is the gas den- 
sity at Ri. We assume that the density distribution is 
fixed a priori with a disk thickness parameter H oc R and 
H/R = 0.12. For HR 4796A, we further assume R t = 
60 AU and an outer disk radius R = 80 AU. The disk 
thickness parameter H = 8.4 at 70 AU is approximately 
consistent with our estimate for an atomic gas scale height, 
H = 8.1 AU (see Section 4). We fix the dust mass M^ust 
= 0.25 M® and calculate the gas chemistry and temper- 
ature distribution assuming gas:dust ratios of 2, 10, and 
100. An additional gas:dust ratio = 100 model is calcu- 
lated assuming a constant escape probability, (3=1. The 
model parameters are summarized in Table 4. 

The calculated temperature structures for the HR 
4796A circumstellar disk models are shown in Figure 3. 
Deep inside the disk, [O I] and [C II] cooling balance pho- 
toelectric heating and heating by H 2 formation. With in- 
creasing gas mass, fine-structure line cooling becomes more 
efficient while photoelectric heating remains constant. The 
photoelectric heating rate does not change because the 
dust mass is fixed in all of the models. Thus, more massive 
disk models possess significantly cooler top layers. Since 
the one dimensional escape probability is only a crude ap- 
proximation to the actual escape probabilty, we also model 
the HR 4796A disk with (3 = 1 for a gas:dust ratio of 100 
(in which the [O I] 63 /zm line has an optical depth of ~10 
at 80 AU). In addition to escape along the line of sight to 
the star, line radiation can also escape perpendicular to 
the disk and through the outer radius. Therefore, mod- 
els which include the one dimensional escape probability 
approximation yield an upper limit to the disk midplane 
gas temperature while models with (3=1 yield a lower 
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limit. The gas temperatures estimated with and without 
the one dimensional escape probability approximation dif- 
fer by less than a factor of 2 or by less than 30% in the 
disk thickness parameter. 

In Table 5, we list the predicted gas temperatures aver- 
aged over one scaleheight at 70 AU and the predicted H I, 
H 2 , Zn II, C II and CO gas column densities assuming a 
linc-of-sight inclined 17° from the disk midplane. The col- 
umn densities estimated from these detailed heating and 
cooling models suggest that the measured column densities 
are not constraining, that the gas:dust ratio = 100 model 
is allowed for all species, except for C II. The measured 
C II column density, N(C II) = 1.4 x 10 14 cm~ 2 , is most 
consistent with the M gas = 0.5 M® model calculated using 
the 1-D escape probability approximation which predicts 
N(C II) = 1.5 x 10 14 cm- 2 . This result is consistent with 
the simple analytical estimate for the gas mass, M gas < 
1.0 M® 

6. DISCUSSION 

In general, the simple analytic model and the detailed 
chemical model estimate similar temperatures and gas 
masses for the circumstellar disk around HR 4796A. The 
discrepancy between the predicted Zn II column densities 
from the simple analytical model and the detailed chemi- 
cal model may be due to (1) the lack of detailed structure 
in the simple model and (2) the different surface density 
distributions used. In the simple model, the entire disk is 
assumed to be either atomic or molecular even though the 
disk possesses a more complicated chemical structure. For 
the low gas:dust ratio case, the upper regions of the disk 
are predominantly atomic while the midplane consists of 
cooler molecular material. The line-of-sight intercepts the 
upper layer of atomic hydrogen but few molecules of cooler 
H2. The surface density, S, in the simple model is a Gaus- 
sian superimposed on a R~ 2 5 power law while the surface 
density in the chemical model is only described by a R~ 2 5 
power law. Despite the difference in assumed functions, 
the surface densities in the models agree to within a factor 
of two at all radii in the disk suggesting that this difference 
can not account for the discrepancy in the predicted Zn II 
column densities. 

The circumstellar disk around HR 4796 probably pos- 
sesses <1.0 M® H I which corresponds to a gasrdust ratio 
<4.0 although a gas:dust ratio of 100 can not be defini- 
tively ruled out. The atmosphere of a giant planet con- 
tains approximately 275 M® of H 2 (Tholen, Tejfel, & Cox 
2000). This is substantially more than either <1.0 M® cir- 
cumstellar gas we infer to exist around HR 4796A or even 
the <5.0 M® gas that would exist if the gas:dust ratio were 
interstellar. Therefore, there is currently not enough cir- 
cumstellar gas around HR 4796A to form the atmosphere 
of a Jovian planet. If Jupiters exist around HR 4796A, 



they must have formed within ^8 Myr, significantly less 
time than is required by multi-stage accumulation (10 - 
100 Myr). 

7. CONCLUSIONS 

We have obtained ultraviolet spectra (905 A - 1187 A, 
1425 A - 1620 A, 1880 A - 2155 A) of the nearby, debris 
disk object HR 4796A using FUSE and STIS on HST. We 
argue the following: 

1. The X-ray flux from the M-dwarf companion HR 
4796B does not significantly heat circumstellar gas around 
HR 4796A. The estimated X-ray heating rate is at least 
<~40 times smaller than the estimated photoelectric heat- 
ing rate. 

2. Circumstellar gas around HR 4796 A may be warmed 
to a temperature T gas <~ 65 K by photoelectric heating 
if the gas:dust ratio is 4, suggesting a molecular gas disk 
thickness parameter H — 6.1 AU and an atomic gas disk 
thickness parameter H — 8.1 AU at 70 AU from the star. 
This is approximately consistent with detailed heating and 
cooling models which suggest a maximum gas tempera- 
ture, T gas 61 K, for a gas:dust ratio of 2, assuming a 1-D 
approximation for the escape probability of line photons. 

3. The circumstellar disk around HR 4796 A possesses 
very little H 2 . We estimate that <1.1 x 10~ 3 M® H 2 ex- 
ists around HR 4796A assuming that gas is in hydrostatic 
equilibrium, is well-mixed, and has a temperature T gas = 
65 K. 

4. We estimate that <1.0 M® atomic hydrogen exists 
around HR 4796A. This corresponds to a gas:dust mass ra- 
tio <4.0 assuming that gas is in hydrostatic equilibrium, 
is well- mixed, and has a temperature T gas = 65 K. 

5. Measurements of C II column density may provide 
the strongest constraint on the gas mass in the HR 4796A 
circumstellar disk. While the measured H 2 , CO, and Zn II 
column densities are consistent with the column densities 
predicted for the gas:dust ratio = 100 model, the measured 
C II column denisty N(C II) = 1.4 x 10 14 cm~ 2 is most 
consistent with the gas:dust ratio = 2 model. 

6. The circumstellar disk around HR 4796 A is too 
gas depleted to support the formation of giant planets. 
Any Jovian planets in this system must have formed on a 
timescale <8 Myr. 
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ministration (NASA) and was partially carried out at 
the Jet Propulsion Laboratory which is managed by the 
California Institute of Technology under a contract with 
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Table 1 
HR 4796A Properties 

Quantity Adopted Value Reference 



Primary Spectral Type 


A0V 


1 


Distance 


67.1±3.4 pc 


2 


Effective Temperature {T e ff) 


10,000±500 K 




Stellar Radius 


1.7 R Q 


2 


Stellar Luminosity 


21 L & 


2 


Stellar Mass (mJ) 


2.4 M Q 




Rotational Velocity (vsini) 


152±8 km/sec 


3 


Fractional Dust Luminosity 






(Lm/L.) 


5xlCT 3 


2 


Estimated Age 


8±2 Myr 


4 


Grain Temperature 


110 K 


2 


Dust Distance (i?o) 


70±1 AU 


5 


Dust Mass (M dust ) 


>0.25 M e 


6 



References. — (1) Hoffleit & Warren (1991); (2) Jura et al. 
(1998); (3) Royer et al. (2002); (4) Stauffer et al. (1995); (5) 
Schneider et al. (1999); (6) Greaves et al. (2000) 
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Table 2 
Interstellar Gas Properties 





Ground State 






Species Wavelength 


Energy W\ 


N 


V 


(A) 


(cm -1 ) (mA) 


(cm -2 ) 


(km/sec) 



Si II 


1526.14 


0.00 


65 


>2.4xl0 13 


-5.0 


Si II 


1526.14 


0.00 


20 


>7.4xl0 12 


-14.4 


Fe II 


1608.45 


0.00 


14 


>9.9xl0 12 


-5.5 


Fe II 


1608.45 


0.00 


15 


M.lxlO 13 


-14.6 


Zn II 


2026.14 


0.00 


7 


>3.8xlO n 


-6.0 



Table 3 

clrcumstellar gas upper limits 



Ground State Gas:Dust 
Species Wavelength Energy W\ N M gas M H Ratio 

(A) (cm -1 ') (mA) (cm- 2 ) (M e ) (M e ) 



CII 




1036.33 


0.00 


<160 


<1.4xl0 14 


<2.4xl0" 


-5 


<4.5xl0~ 


-3 


<0.018 


OI 




1039.23 


0.00 


<250 


<2.8xl0 15 


<6.4xl0" 


-4 


<0.051 




<0.20 


Zn II 




2026.14 


0.00 


<5 


<2.6xl0 12 


<2.4xl0" 


-6 


<1.0 




<4.0 


H 2 (J= 


=0) 


1036.54 


0.00 


<260 


<1.0xl0 15 


<2.5xl0" 


-4 


<2.5xl0" 


-4 


<1.0xl0" 3 


H 2 (J= 


=1) 


1038.16 


118.16 


<300 


<3.7xl0 15 


<9.1xl0" 


-4 


O.lxlO" 


-4 


<3.6xl0- 3 


CO (J= 


=1) 


1544.45 


3.85 


<4 


<1.2xl0 13 


<4.1xl0" 


-5 


N/A 




N/A 


CO (J= 


=2) 


1544.39 


11.58 


<4 


<2.4xl0 13 


<8.3xl0" 


-8 


N/A 




N/A 



Table 4 
HR 4796A Model Parameters 



Model 


M gas 


Gas:Dust 


a 


escape approach 




(M ffi ) 


Ratio 


(cm 2 /H-atom) 




Model 1 


0.5 


2 


1.17 10- 21 


TH85 


Model 2 


2.5 


10 


2.34 10- 22 


TH85 


Model 3 


25.0 


100 


2.34 10- 23 


TH85 


Model 4 


25.0 


100 


2.34 1Q- 23 


(3 = 1 



Table 5 

Model Predicted Gas Temperatures and Column Densities 

Averaged Predicted Predicted Predicted Predicted Predicted 
Model T ffas (70AU) A H A H2 A z „+ JV C + A co 

(K) (cm- 2 ) (cm- 2 ) (cm -2 ) (cm -2 ) (cm -2 ) 



Model 1 61 1.1 x 10 18 2.3 x 10 13 4.2 x 10 10 1.5 x 10 14 4.2 x 10 6 

Model 2 45 5.3 x 10 18 8.4 x 10 13 2.0 x 10 11 7.5 x 10 14 1.1 x 10 8 

Model 3 52 5.3 x 10 19 1.0 x 10 15 2.0 x 10 12 7.3 x 10 15 9.9 x 10 9 

Model 4 22 
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Fig. 1 — The interstellar (a) Si II A 1526.71 and (b) Fe II A 1608.45, and (c) Zn II A 2026.14 lines obtained with STIS. The vertical lines 
indicate the wavelengths of the Si II, Fe II, and Zn II lines at the star's velocity.. 
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20 40 60 80 100 

Temperature (K) 

Fig. 2. — The gas temperature, estimated from equation (3) plotted as a function of the gas:dust ratio or gas mass (dotted line). The 
circumstellar gas mass upper limit inferred from the Zn II line at 2026.14 A plotted as a function of gas temperature (solid line). In these 
simple models, the gas has a constant temperature throughout the disk. Both relationships constrain the circumstellar gas mass as a function 
of temperature. 
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M(dust) = 0.25, M(gas) = 0.5 M(dust) = 0.25, M(gas) = 2.5 




r [AU] r [AU] 

Fig. 3. — (a)-(c) The gas temperature structure in a slice perpendicular to the disk for models with gasidust ratios of 2, 10, and 100, 
calculated using the one dimensional escape probability formalism discussed in the text, and (d) the temperature structure for a model with 
a gas:dust ratio of 100 calculated using f3 = 1.0 



